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Abstract. We study mobility and interaction of gap solitons in a Bose-Einstein 
condensate (BEC) confined by an optical lattice potential. Such localized wavepackets 
can exist only in the gaps of the matter-wave band-gap spectrum and their interaction 
properties are shown to serve as a measure of discreteness imposed onto a BEC by 
the lattice potential. We show that inelastic collisions of two weakly localized near- 
the-band-edge gap solitons provide simple and effective means for generating strongly 
localized in-gap solitons through soliton fusion. 
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1. Introduction 

Periodic potentials created by optical lattices are by now recognized as a powerful tool 
for controlling and manipulating nonlinear matter waves through dynamical diffraction 
management plEj- One of the important implications of such management is creation 
of nonlinearly localized matter-wave packets, the so-called bright atomic solitons, in 
a Bose-Einstein condensate with a positive scattering length, i.e. with repulsive inter- 
atomic interaction. Such a counter-intuitive effect can occur as wave localization in the 
gaps of the band-gap spectrum due to the Bragg scattering of coherent matter waves by 
the periodic potential of an optical lattice Since the repulsive atomic interactions 
prohibit a collapse of the bright wavepackets and the solitonic nature of their evolution 
ensures fixed spatial phase variations due to mean-field effects, the atomic gap solitons 
may represent an attractive high-density source for atomic interferometry 

The possibility of nonlinear localization in a repulsive BEG was established 
theoretically IH El El E!, being recently confirmed experimentally in the case of one- 
dimensional optical lattices [S]. The experimental challenge of the gap soliton generation 
and manipulation is two- fold. First, the BEG wavepacket, initially loaded into a ground 
(Bloch) state of the optical lattice potential, should be accelerated to the edge of the 
first Brillouin zone. The regime of the wavepacket preparation and the dynamical 
transition to the band edge dramatically affects the outcome of the experiment 
and in the best case scenario only weakly localized low-atom-number gap solitons near 
the bottom edge of the spectral gap can be generated using the current experimental 
techniques |H]. The depth of the gap, where the gap solitons contain large atom numbers 
and are well localized is yet to be accessed. Secondly, the possibility of the gap 
soliton manipulation, which is essential for interferometric applications, can be limited 
by impaired mobility of the gap solitons due to effects of lattice discreteness |TT] . 

In this paper, we examine the problem of mobility and interactions of gap solitons 
in a one-dimensional optical lattice without an additional harmonic confinement, within 
the framework of the continuous Gross-Pitaevskii model with a periodic potential. We 
show that, as the chemical potential of a BEG scans the spectral gap, the gap solitons 
display a variety of scattering properties ranging from those of "ideal" solitons without 
a periodic potential to those of nonlinear localized excitations of discrete nonlinear 
lattices [12 . In addition, we suggest that the inelastic collisions of near-the-band- 
edge gap solitons can provide reliable means for generating highly-localized high-density 
immobile atomic solitons in the depth of the spectral gap. 

2. Model 

We consider the dynamics of a cigar-shape BEG cloud in the presence of a one- 
dimensional optical lattice. In dimensionless units the Gross-Pitaevskii (GP) equation 
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governing the evolution of the condensate wavefunction can be written as: 



K4l^^^<^' + ^'"0"• 



where the optical lattice potential V{x) = Vosiri^ {nx / d) is uniform and characterized 
by its depth Vq measured in units of lattice recoil energy Er = h'^n'^ /{2md'^), and lattice 
period d. The condensate wave function is normalized to: dx\ip{x,t)\'^ = Nqid, 
where is the number of atoms, gm = 2|a<j|/ao, a<j is the s-wave scattering length, and 
Oq the characteristic harmonic oscillator length in the direction of tight confinement. 
For ^''Rb atoms with m = 1.44 x 10~^^kg, = 5.3nm and assuming a transverse 
frequency u± ^ 271 x lO^Hz, gives gin ~ 10"^. The coefficient a = sgn{as) = ±1 
characterizes the type of the atomic interactions. In what follows we consider repulsive 
interaction, i.e. cr = 1, and set d = 1, which corresponds to ~ 1.1 fim lattice spacing. 

Stationary states of a condensate are described by solutions of Eq. (jT)) of the form: 
ilj{x;t) = (f){x) exp{—ifit), where the steady-state wave function (j){x) obeys the time- 
independent GP equation 



In the case of a noninteracting condensate {a = 0) stationary solutions of equation (j2I) 
are Bloch waves, ~ b{x, k)exp{ikx), where b{x) have periodicity of the lattice 

and k is the quasi-momentum. The linear spectrum n{k) of the Bloch states has 
a characteristic band-gap structure, with the lowest two bands shown in figure Ufa) 
in the parameter plane (Vq, /i). The shaded areas correspond to the regions where 
oscillating solutions of equation Q (Bloch waves) exist. The clear areas correspond 
to the spectral gaps which appear due to the Bragg scattering of the matter waves in 
a periodic structure, being the forbidden domains for matter wave propagation. The 
gap below the first band is the trivial semi- infinite "total internal refiection" gap. The 
first finite gap is where the formation of gap solitons in repulsive BEG occurs [HI EI- 
The condition fi = Vq roughly delineates between the regimes of tight-binding (/i < Vq) 
and superfiuid (// > Vq) behavior. In what follows, we investigate the dynamics of gap 
solitons in the lattice of height Vq = 10 which places the first gap into a relatively 
tight-binding regime. 

The nonlinear localization of a BEG with positive scattering length (a = +1) in the 
first spectral gap has been well studied both theoretically [HllZl and experimentally [S]. 
The stationary localized solutions of equation (j2)) - gap solitons - can be found at any 
value of /i inside the spectral gap. The number of atoms contained in a localized state 
is small near the lowest gap edge and raises as the chemical potential moves inside 
the gap inj. The lowest-order families of on-site (OS) and inter-site (IS) gap solitons, 
centered on the lattice minimum or maximum, respectively, are shown in figure ^^b). 
The characteristic "staggered" spatial structure of both OS and IS solitons can be seen 
in figure mc). 




(2) 
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3. Mobility of gap solitons 

In order to study interactions of gap solitons, one has to generate mobile localized states. 
It has been shown, however, that the mobility of a matter wave soliton in an optical 
lattice is restricted by the existence of the so-called Peierls Nabarro (PN) potential bar- 
rier Ej- The concept of the PN barrier originates from the theory of dislocations 
and later studies of mobility of localized nonlinear excitations in discrete lattices [Tni ll4j. 
The PN barrier is introduced as the height of the effective periodic potential generated by 
the lattice discreteness, and it defines the minimum energy required to move the centre 
of mass of a localized wavepacket by one lattice site fHj- Extending this definition to 
the continuous lattice model, one can define the PN barrier as the energy difference 
between the on-site and inter-site localized stationary states |TTj. Since the OS and IS 
states can be mapped to represent a moving mode at different times, the PN barrier can 
be calculated as the smallest amount of energy that a gap soliton needs to gain in order 
to start moving along the lattice. Even in the tight-binding regime, the amount of the 
energy needed to initiate the soliton motion is not in agreement with the estimate of the 
PN barrier height derived from the corresponding discrete nonlinear Schrodinger model 

therefore the use of the continuous model in the study of gap soliton motion and 
interactions is well justified. 

A soliton solution of equation Q is characterized by its dynamical invariants: 
number of atoms and energy. The soliton energy can be calculated as: 



The on-site (OS) gap sohton is a localized state with the lowest energy, and therefore it 
is this state that is eventually generated from an arbitrary initial wavepacket. The PN 
potential height associated with this state is given by Upn = Ejs — Eqs- The threshold 
value of the momentum (per atom) that a soliton must have to overcome the barrier 
is defined as kp^ = {2Upn /Ny/'^ . Such a momentum can be applied experimentally 
either by ultrashort laser pulses or a tilt of the optical lattice potential imposed by 
a gravitational field or by a slope of a harmonic trap. The PN potential per atom 
calculated for the family of OS gap solitons within the first gap is shown in figure Efa). 
In contrast, the inter-site state corresponds to a maximum of the PN potential and is 
unstable against small perturbations of its center-of-mass position. Such perturbation 
will lead to conversion of the IS soliton into an OS soliton, as shown in figure El 

To investigate the dynamics of the OS gap soliton with a finite velocity, we solve 
the evolution equation numerically. The initial condition is given by ?/'(x; t = 0) = 
(poi^x) exp {—ikox), where 0o is a stationary gap soliton wavefunction, and /cq is an initial 
momentum per atom. The values of were chosen below, above, or at the threshold 
value kpN- The results of the calculations are shown in figure |2tb,c,d). Due to 
the large negative group- velocity dispersion pH], solitons near the edge of the lower 
band exhibit anomalous steering, i.e. initially the direction of motion is opposite 
to the momentum kick received. Oscillating motion within a single lattice site has 




(3) 
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been observed for an initial ko values below and at the threshold [see figure |2fb,c)]. 
Above the threshold kinetic energy (which is still an order of magnitude lower than 
an individual lattice well depth) a soliton undergoes free motion across the lattice, as 
shown in figure Efd). During this motion, an atomic soliton evolves through many states 
centered on different points of the optical lattice. Periodic revival of the states closely 
resembling the stationary states of the IS and OS families [figure HJ^c)] can be observed, 
however due to a loss of atoms neither number of atoms contained in a soliton nor /i 
are conserved during the evolution. The loss of atoms from the moving localized state 
increases with growing ko. 

The mobility properties of gap solitons vary depending on their chemical potential. 
Close to the lowest gap edge (yU < 7.5 in our case), gap solitons are well described by 
the effective nonlinear Schrodinger equation (NLS) for the slowly varying envelope of a 
Bloch state [H]. In this regime, solitons are expected to behave like lattice-free envelope 
solitons of the integrable NLS model. Indeed, they exhibit free motion for arbitrary 
small initial momentum due to very small energy difference between the OS and IS 
stationary solutions of equation (0). These states are weakly localized and contain only 
about 100 atoms. Deeper in the gap (7.5 < ji < 8.2), the mobility properties of solitons 
are similar to those of localized solutions of the discrete NLS equation, and can be very 
well described by the concept of the PN barrier. In the regime where the solitons are 
strongly localized in the vicinity of a single lattice well (yU > 8.2), their dynamics is 
dominated by oscillations of the centre of mass within a single lattice site for moderate 
initial perturbations. The soliton motion with a large initial velocity is accompanied by 
significant loss of atoms from the localized structure. 

4. Binary collisions of gap solitons 

Depending on the value of the chemical potential and degree of localization, gap solitons 
display a variety of scattering properties. In particular, we study binary collisions 
of in-phase solitons numerically for a range of values of the chemical potentials near 
the bottom edge of the gap, which is the experimentally accessible region of soliton 
generation (S]. We keep the initial separation of soliton centers the same in all cases and 
chose the same initial velocity above the PN threshold. Similar interaction conditions 
can be realized by imposing an additional harmonic confinement onto the system. 

Three regimes of soliton collisions have been identified in our simulations. The first 
one is the regime where in-phase solitons (A(/) = 0) interact elastically, i.e. without 
momentum and energy exchange [see figure |3 (a)]. In this narrow parameter domain 
(7.3 < fi < 7.35) the dynamics of the gap solitons resembles the dynamics of "ideal" 
solitons of the completely integrable NLS equation without the lattice potential. The 
total momentum of the system is defined as: 




(4) 
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and its zero value before and after collision indicates that the soliton collisions are not 
affected by the lattice periodicity. 

Deeper in the gap the solitons start to experience stronger localization and their 
dynamics become inelastic. Significant radiation usually accompanies the collisions, 
and the total momentum of the system is no longer conserved. This is the typical 
signature of a discrete system [T2j, and a consequence of nonintegrability introduced by 
a periodic potential. The effect of small perturbations (numerical noise which emulates 
noise in the real physical system) accounts for the spontaneous symmetry breaking (SB) 
and generation of finite momentum during the collision process ^2)- In this "weakly 
discrete" regime of the soliton dynamics, the solitons can exhibit multiple collisions |12j . 
In the multiple collisions solitons "bounce" from each other leading to either final 
separation or fusion, see figure |ll^b,c). For /i > 7.43 the solitons merge into a single 
localized state, as seen in figure Efd), which due to SB moves in an unpredictable 
direction across the lattice and exhibits large amplitude oscillations. The appearance 
of SB, multiple collisions and the total momentum generation has been observed for 
a range of the chemical potential 7.35 < /i < 7.5. 

For our choice of the optical lattice depth, the solitons relatively close to the 
lower gap edge (e.g. for /i = 7.5) already exhibit strongly inelastic scattering behavior 
characteristic of a discrete system. For the values of the chemical potential ^ > 7.5 
the binary collisions lead to fusion of two solitons into a single localized state. Unlike 
the outcome of fusion after multiple collisions, the resulting solitons move with nearly 
zero velocity and damped amplitude oscillations, as can be seen in figure El Location 
of the initial states (containing around 170 atoms each) on the branches of stationary 
solutions, is indicated by point A in figure ^b). After collision, an immobile highly- 
localized soliton containing a larger number of atoms is generated. From the vs. /x 
dependence [see figure Hfb)] it can be inferred that the generated soliton, consisting 
of approximately 270 atoms, lies deeper in the gap than the initial states [see point 
B in figure ^b)]. Therefore, inelastic collision of weakly localized near-the-gap-edge 
solitons can, in principle, be employed for a formation of fundamental atomic gap 
solitons with a larger number of atoms in the depth of the spectral gap, i.e. in the 
region of parameters that is currently inaccessible experimentally. This technique can 
be useful for any number of colliding solitons near the band edge. In particular, the case 
of triple-soliton collisions is illustrated in figure were the phase of the central soliton 
was inverted by vr to achieve soliton merging. 

For solitons located even deeper in the gap, the binary collisions inevitably result 
in their fusion, and the emerging localized states remain confined within a single lattice 
well, although their centre of mass can oscillate around the potential minimum. The 
amplitude of the centre of mass oscillations increases with /i, however the value of the 
total momentum averaged over the time At at ~ 20 oscillations: {p) = J^°~^^^ pdt, 
remains close to zero. The values of the post-collision velocity of the localized state 
(relative to the initial velocity vq) and the average momentum vs. chemical potential 
of the colliding solitons in the vicinity of the gap edge are plotted in figure H It can 
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be clearly seen that the initial stage of quasi-elastic collisions {{p) ^ 0, v = vq) and 
final stage of the soliton fusion ((p) «i 0, f ~ 0) is separated by the symmetry breaking 
stage, where the outcome of collisions is not predictable. We note that the symmetry 
breaking outcome of soliton collisions is also characteristic of binary collisions of solitons 
with non-zero phase difference, as in the case of a localized state in the discrete NLS 
model [12] . 



5. Conclusions 



We have studied numerically both mobilility and interaction of the matter-wave gap 
solitons in an experimentally accessible regime, i.e. near the lowest edge of the spectral 
gap of an optical lattice. For the collision velocities above the Peierls-Nabarro mobility 
threshold, the gap solitons demonstrate scattering dynamics scenarios ranging from 
quasi-elastic collisions to soliton fusion. By employing the inelastic collision of two or 
three gap solitons in configurations that can be realized experimentally, for instance, in a 
combination of an optical lattice and a harmonic trap, we showed that the soliton fusion 
can be successfully used to generate stationary matter-wave solitons deep in the spectral 
gap. Applying the same technique to solitons with different amplitudes can potentially 
enable generation and steering of a mobile high-density soliton. Such a process of gap- 
soliton manipulation may eventually lead to new possibilities for the preparation of 
sources for atomic interferometry 
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Figure 1. (a) Band-gap structure of the matter-wave spectrum /i(Vb). (b) Families of 
the on-site (OS, dotted) and inter-site (IS, solid) gap solitons for Vq = 10.0, within the 
first gap. Point A (B) corresponds to the location of an initial (final) state in figure\^ 
(c) Spatial profiles of the OS (top) and IS (bottom) gap solitons at = 7.5. 




Figure 2. Top: (a) PN potential for gap solitons presented in figure^^b) as a function 
of fi within the first gap. Bottom: Examples of soliton centres of mass dynamics for 
fi = 7.7 with the PN threshold momentum kp^ = 1.20 x 10^^ and the initial momentum 
(b) below the PN threshold, fco — 6.01 x 10^^, (c) near the threshold ko = 1.20 x 10~^, 
and (d) above the threshold ko = 1.24 x 10~^. 



Interaction of matter-wave gap solitons in optical lattices 



10 




1000 2000 3000 

t 



Figure 3. Evolution of the centre of mass of the initially stationary IS gap soliton at 
fi = 7.7. Insets: (a) Initial spatial profile, and (h) spatial profile at the time t = 3x 10^, 
showing the charactersitic structure of the OS state. 
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Figure 4. Examples of soliton interactions in the nonintegrable NLS regime for 
different values of the chemical potential (a) fi = 7.35; (b) /i = 7.405; (c) fi ^ 7.41; 
(d) n = 7.45 and ko = 0.1. 
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Figure 5. Top: Typical scenario of soliton fusion with zero final velocity in the discrete 
NLS regime (fj, = 7.5) Initial momentum of the soliton pair |fco| — 0.1.; Inset: Solid - 
typical spatial profile of the post- collision trapped state at t = 10'^ , dashed - exact OS 
solution corresponding to point B in figure^fb), fi — 8.36. Below: Evolution of the 
peak density of the trapped state showing transition to a trapped highly-localized in-gap 
state corresponding to the marked point B in Roure\lVb). 
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Figure 6. Typical scenario of triple- soliton fusion in the discrete NLS regime 
(j-i ~ 7.5j, with the central immobile soliton -K-out of phase with the others. Initial 
momentum of the off-centre soliton pair jfepl ~ O.2.; Inset: typical spatial profile of the 
post-collision trapped state at t = 10^. 
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Figure 7. Post-collision (at t = 10^) (a) time- averaged momentum of gap solitons 
as a function of jj. within the first gap and (b) relative velocity of the soliton pairs. 
Behavior characteristic for fco = 0.1. 



